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ABSTRACT 

We construct a map of the time derivative of the gravitational potential traced by SDSS Luminous 
Red Galaxies. The potential decays on large scales due to cosmic acceleration, leaving an imprint 
on cosmic microwave background (CMB) radiation through the integrated Sachs- Wolfe (ISW) effect. 
With a template fit, we directly measure this signature on the CMB at a 2-cr confidence level. The 
measurement is consistent with the cross-correlation statistic, strengthening the claim that dark energy 
is indeed the cause of the correlation. This new approach potentially simplifies the cosmological 
interpretation. Our constructed linear ISW map shows no evidence for degree-scale cold and hot 
spots associated with supervoid and supercluster structures. This suggests that the linear ISW effect 
in a concordance ACDM cosmology is insufficient to explain the strong CMB imprints from these 
structures that we previously reported. 

Subject headings: cosmic microwave background — cosmology: observations — large-scale structure 
of universe — methods: statistical 



1. INTRODUCTION 

Large-scale structures in the low-redshift Universe 
leave a mark on the cosmic microwave background 
(CMB) radiation th rough the late-time i ntegrated Sachs- 
Wolfe (ISW) effect dSachsfc Wolfe|1967| ). In ACDM, the 
expansion of the Universe accelerates at late times, caus- 
ing gravitational potentials to decay. The effect gives 
an energy boost to photons traveling through massive 
structures and degrades the energy of photons crossing 
under-dense voids. In a flat universe, the effect occurs 



only in the presence of dark energy ( Crittenden & Turok 
1996). The evolution of the ISW signal provides a con- 
stramt on the dark energy equation of state, and is a 
unique probe of cosmology, independent from studies us- 
ing supernovae as standard candles, the bary on acoustic 
scale, galaxy cluster counts, or weak lensinj 
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et al.| |2005| _ 

The linear ISW signal on the CMB is difficult to mea- 
sure because the temperature fluctuations from the ISW 
effect are about an order of magnitude smaller than those 
of the primary CMB. However, the effect can be mea- 
sured statistically by correlating the structure at low red- 
shift with the CMB temperature. This has been carried 
out with many galaxy sa mples with redshift z < 2, in - 
cluding the APM surve y (Fosalba & Gaztahaga 2004), 
2MASS (lAfshordi et al. 12004; Rassat et al.||2007D, SDSS 



Luminous Red Galaxies (Scranto n et al. | [ 2003 " [Fosalba 
let al 120031 |Padmanabhan et al. | 2005p , SD SgMain Sam- 
ple (|Cabre et al.Jj2006[ ), SPSS quasars (|Giannantonio 



eT!d1|2006l), NVS 
2004} |lNolta et al 



2004 Raccanel 



radio sources (Boughn & Critt enden] 
et al. 2008) and the 



X-ray background ( jBoughn & Cr ittenden 2004). The 
correlation has been studied using wavelet analyses as 
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well ( McEwe n et al.|2007| |Pietrobon et al.|2006| [McEwen 
et al. 2008| ). Most recently, the combined analys is of mul- 



tiple surveys gives a detection s ignificance > 4a (Ho et al. 



2008) |Giannantonio et al.|2008| >. 

On non- linear scales, gravitational potentials evolve 
through structure formati on processes, giving rise to the 
Rees-Sciama (RS) effect ([Rees fc Sciama 1968 |Sakai fc| 



Inoue|[2008j). The RS effect dominates over the linear 
ISW signal a t small angular scales (multip oles > 1000) 
on the CMB ( |Cooray|2002bl |Cai et al.p009| > . Other sec- 
ondary anisotropics on the CMB include the Sunyaev- 
Zeldovich (SZ) effect induced by inverse- Compton scat- 
terin g of CMB photons by hot gas in massive clus- 
ters (Sunyaev & Zeldovich 1972). Gravitational lens- 
ing by foreground structures also induces a temperature 
anisotropy, although both of th ese effects are re stricted to 
arcminute scales on the CMB ( |CoorayJ|2002aj ). On high 
redshift galaxy samples, the interpretation of degree scale 
matter-CMB correlations can be complicated b y magni- 
fication bias, which boosts the ISW signal ( |Loverde et al.| 
2007), a s well as large-scale correlations ar ising at reion- 
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ization ( Giannantonio & Crittenden 

In this paper, we construct an ISW map from a sam- 
ple of luminous red galaxies (LRGs), assuming linear 
growth of density fluctuations, and using a parameter- 
free Voronoi tessellation technique to add the potential 
from each galaxy directly. We measure the ISW signal 
on the CMB using the constructed map as a template. 
Typically, the ISW measurement is degraded by cosmic 
variance on the CMB temperature as well as variance in 
the local large-scale structure. Our approach utilizes in- 
formation about the observed density field traced by the 
LRG survey to remove the effects of local cosmic vari- 
ance, potentially improving the detection. We compare 
this technique with the cross-correlation statistic. For 
the LRG sample, the two measurements are consistent, 
although the error analysis is simplified considerably by 
the template fit approach because a Monte-Carlo analy- 
sis is not required. 

Template fits (or matched filter analyses) are com- 
monly employed in studying Galactic foreground emis- 
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sion on the CMB (de Oliveira-Costa et al. 1999), and 
have been u sed in the context of the SZ effect ( |Hansen 



et al. 2005). The application to ISW ha s been pro 
posed by a numbe r of groups, including [He rnandez- 
Monteagudo| d2008| ), iFYommert et aTJ ([2008| , an d [Bar 
reiro et al.| (|2008p. I Hernandez- Mont eagudol and prom- 



IT? 

mert et aTf find that by removing cosmic variance from 



the measurement, the ISW detection significance can im- 
prove by 10% over the correlati on function measure ment, 
a conclusion also supported by Cabre et al. ( 2007| ). 

A map of the foreground anisotropics is also of interest 
in itself. Large structures, especially at low redshift, can 
potentially produce significant anisotropi cs on the CMB 
( |Maturi et aT]|2007| [Sakai fc Inoue||2008| ). These might 
explain CM B anomalies includi ng the alignment of low 
multipoles fllnoue fc S ilk 200 7[), and the 5° Cold Spot 
flRudnick et al.||2007| ). |Granett , etH] ( |2008D found a 4- 
cr-significant signature corresponding to supervoids and 
superclusters in in the Wilkinson Microwave Anisotropy 
Probe (WMA P) maps, which hin ts at physics outside of 
ACDM (e.g. |Hunt fc Sarkar|2008[ ). We will further inves- 
tigate this signal again using template fitting techniques. 

Unless noted, we employ joint WMAP 5, super novae 
and BAO ACDM paramete rs: fl c = 0.233, Q b = 0.0462, 
h = 0.701 with cr 8 = 0.817 ( |Komatsu et~aL||2009[ ). 



2. DATA 
2.1. Luminous Red Galaxies 

We b ase our study on the Sloan Dig ital Sky Survey 
(SDSS) flAdelman-McCarthy et~aL|2008| ) Luminous Red 
Galaxy (LRG) sample. The photometric galaxy catalog 
traces large-scale structure to z = 0.8 over 7500 square 
degrees of contiguous sky. LRGs are elliptical galaxies in 
massi ve galaxy clusters representing large dark- matter 
halos ( |Blake et al.||2008] ), and are thought to b e phys- 
ically similar obj ects across their redshift range ( |Eisen-| 
|stein et al.||200l| |Wake et aLl|2006| ). This makes them 
excellent, albeit sparse, tracers of the cosmic matter dis- 
tribution on scales > 10 Mpc. Our sample is designed 
to match the spectroscopic targets in the 2SLAQ LRG 
survey ( |Cannon et al.||2006 ). In particular, we apply the 
stringent i < 19.8 and d perp > 0.55 cuts. Within the 
DR6 imaging area about the North Galactic Pole, the 
catalog includes 746,962 objects; for our fiducial sample, 
though, we use the 400,000 closest to the median in red- 
shift. Phot ometric redshifts are obtained from |Oyaizu| 
et al.|(2008) we use the results from the CC2 algorithm. 
The median photometric error of the sample is <j z = .039. 
The redshift distribution is plotted in Fig. [I] It extends 
from 0.45 < z < 0.75 with median z = 0.52. A map of 
the projected density is shown in Fig. [3J 

We find a systematic trend in the photometric redshift 
distribution with Galactic latitude; see Fig. [I] At low 
latitude, the distribution shifts toward higher redshift. 
Although the shift in the median redshift is smaller than 
the typical redshift error, it produces a strong gradient 
in the galaxy density at high and low redshift where the 
selection function is steep, which contaminates the 3D 
density reconstruction. We attribute this to residual er- 
rors in the Galactic extinction correction. The trend of 
median redshift with extinction, E(B-V), is shown in Fig. 
|2j Due to this issue, the useable redshift range is limited 
to the peak of the distribution, as we discuss in §4.3. 
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Fig. 1. — The redshift distribution of the LRG catalog binned 
with Az = .01 is plotted. We constructed the photometric catalog 
from SDSS DR6 using the same selection criteria as the 2SLAQ 
spectroscopic LRG survey. The renormalized spectroscopic redshift 
distribution is plotted in the top frame (shaded). Also plotted are 
the photometric redshift distributions of a high Galactic latitude 
sample with b > 68° and a low latitude cut with 20° < b < 40°. 
The two samples contain equal numbers of galaxies. The difference 
in number density between the full and subsamples are plotted in 
the bottom frame. The shaded range represents Poisson errors. 
We find that the low latitude sample has a systematically higher 
median redshift, but is within the redshift error distribution. 
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Fig. 2. — We attribute the shift of the photometric redshift 
distribution with Galactic latitude to residual errors in the extinc- 
tion correction. Shown is the median redshift in 12° bins (Healpix 
pix els at n^^-lG) pl otted against the Galactic extinction, E(B- 
V)( |Schlegel et al.|1998| >. Over-plotted are the mean and one stan- 
dard deviation limits, as well as a linear fit. The correlation co- 
efficient is r = 0.38. This trend is an order of magnitude below 
the size of the redshift errors, but systematically degrades the 3D 
density construction; see §4.3. 



2.2. Microwave maps 

We use WMAP five-year maps, quoting results from 
the the Q, V and W frequency band maps, in which 
the individual differencing assembly maps have been co- 
added. To estimate Galactic foreground contamination, 
we use the foreground-reduced Q, V and W maps, the 
Markov-chain Monte Carlo (MCMC) derived tempera- 
ture map and the Internal Linear Combination (ILC ) 
map, all of which are described by Gold et al. (2009). 
The Q, V and W foreground reduced maps are full- 
resolution maps with best-fit foreground templates sub- 
tracted, preserving the noise properties of the original 
maps. The MCMC map was produced using a Monte 
Carlo joint fit of the CMB temperature, polarization and 
foregrounds. The map is smoothed at 1° full- width half- 
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max resolution. The ILC map is foreground-cleaned us- 
ing a minimum- variance fit of the foreground templates. 
It is smoothed to 1° resolution. We present the results 
from the ILC map even though its noise properties are 
not characterized. To each map we apply the KQ75 
galactic foreground and point source mask. We carry 
out all analyses using the Healpix 5 pixeliz ation scheme 
at a resolution of 55'/pixel (n S id e = 64) ( [Gorski et al. 
2005), at finer resolution than the > 1° scales where 
tnehnear ISW effect is thought to operate. The MCMC 
map is shown in Fig. [3j 

3. ISW MAP CONSTRUCTION 

The ISW-induced temperature is an integral over the 
time-derivative of the gravitational potential <£>. In terms 
of the conformal time r, th e temperature fluctu ation 
Tisw along a line of sight, is (Sac hs fc Wolfe|[l967 ) 



ATisw 

^CMB 



/ 



dr 



dr 



(1) 



In linear theory, S evolves with the growth factor D\(z) 
as S(x,z) = J^^=o) $ z — 0)- I n a fl a t? matter domi- 
nated universe D\ = 1/(1 + but in ACDM, it deviates 
from this at low redshift due to accelerated expansion. 
We obtain the time derivative of the potential assuming 
linear theory, i.e. 



d$(x) $(x) d 



dr (1 + z) dr 



[(l + z)D 1 (z)]. 



(2) 



We model the densit y field in the survey using a 
Voronoi tessellation (e.g. jOkabe et al.|2000 van de Wey- 
gaert & Schaap 2009[ ) . Each galaxy occupies a polyhedral 
Voronoi cell ol points closer to that galaxy than to any 
other. The volume Vi of that cell gives a natural estimate 
of the galaxy's overdensity, Si = V{zi)jVi — \, where V(z) 
is the average volume of a galaxy at redshift z. 

We compute the potential as a direct sum over these 
polyhedral cells of volume Vi and overdensity Si, concen- 
trated at the comoving positions of their galaxies. In 
the Newtonian limit, the potential is 



#(r): 



Gp cr oflm /-, . x ^iVi 

— (l + ^ med ) ^ — - 



m 



(1 + Z me d 



V(zi) - Vi 



(3) 



(4) 



Here, p cr ^ is the critical density at z = 0. In evaluating 
this expression, we put the whole survey at the median 
redshift of the sample, z me d = 0.52, multiplying by a fac- 
tor ^ m ed+i)gi(^med) ^ After computing the potential, we 

then divide & by this factor to get the proper amplitude. 
Explicitly, the observed density is the true density con- 
volved with the photometric redshift error distribution, 
and attenuated by the top-hat survey window function: 

5=(5*f z )xW(r). 

Our method is relatively slow computationally, since it 
involves a sum over all galaxies for each point where the 
potential is evaluated, and a Voronoi tessellation. For 
our fiducial galaxy sample, the computation takes a few 



hours using several CPU's, while a fast fourier transform 
method could take less than a minute on one CPU. How- 
ever, our method easily allows the true (curved) geom- 
etry of the LRG sample to be included, and integrated 
through to z = 0. Also, the contribution to the potential 
from each galaxy is included in a natural way using the 
parameter-free, highly adaptive Voronoi method. 

We added a buffer of galaxies around the survey, and 
filled survey holes with gala xies sampled at the me an den- 
sity at each redshift, as in |Granett et al.] (2008). Hole 
galaxies (comprising ~ 1/300 of the galaxies) were in- 
cluded in the potential sum, but galaxies neighboring 
any buffer galaxy (i.e. affected by the edge) were not. To 
get the mean galaxy volume function V{z), we spline- 
interpolated between galaxy volumes averaged in bins of 



logarithmic width d\og 



10 ' 



1/200. 



To get the predicted ISW map, we integrated the sig- 
nal (Eq. 0) through lines of sight of Healpix pixels. 
We integrated from z = to 1.5 times the farthest dis- 
tance of a galaxy in the survey, sampling the potential at 
10 (a sufficient number, according to convergence tests) 
equally spaced points in each of three regions: in front 
of, within, and behind the sample. In spots, the nonzero 
potential behind and in front of the survey generated 
~50% of the signal. To lessen the noise from galaxies 
happening to lie too close to a point r where the poten- 
tial was evaluated, we softened the potential by setting 
the distance to a galaxy to max(|r — |, 10 ft^Mpc). 
The resulting map is shown in Fig. [3] 

A constant multiplicative bias factor b Q is thought to 
describe the ratio of the galaxy to matter overdensity 
adequately on the large scales of the linear ISW effect: 
S dm — (l/b a )S a . In the case of the LRG sample, b g ~ 2.2 



(cr 8 = 0.82) flBlake et al.p )Q8). We expect the true ISW 
fluctuations ^isw to be lower than our construction by 
a factor of b g : T^^^ = (l/b g )T^^ G \ This dependence 
on the bias allows for b g to be fit directly, independent 
of the power spectrum amplitude. 



4. STATISTICAL SIGNIFICANCE 

4.1. Amplitude fit 

To judge the reality of the reconstructed ISW signal, 
we measure its contribution to the microwave sky with a 
template fit approach. The observed microwave temper- 
ature is the sum of the primary CMB anisotropy, galac- 
tic and extragalactic foreground contributions and the 

ISW signal, T^^ G ^ /b g . The ISW anisotropy arises from 
the structure sampled by LRG galaxies corrected for the 
galaxy bias. We add the predicted ISW template with a 
scaling parameter A, to be fit: 



T = TpHm + Tf g + — Tisw- 

On 



(5) 



Fixing the ISW and foreground anisotropics, the like- 
lihood function, in terms of the pixel CMB covariance 

matrix, Cij — (^prim,i^prim,j)? is, 



A , 



[ ISW. 
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Fig. 3. — The maps used in our analysis. The left frame is the LRG galaxy overdensity within the SDSS DR6 survey footprint. At 
center is the reconstructed ISW signal in /iK. Right is the foreground-cleaned CMB temperature from the WMAP5 MCMC analysis with 
the galactic foreground and point source mask applied. The maps are in celestial coordinates with 55'Zpixel. The grid spacing is 30° with 
RA=180° at the center, and RA increasing to the left. 



and maximizing this over A gives, 



X = b ( 



Y^ij ^ISW,z 



isw j 



The variance of this estimator is, 



i.C-- Ttc 



(7) 



(8) 



This approach can be extended naturally to multiple 
correlated datasets through a joint fit of many template 
maps. The estimator is non-linear in the the ISW tem- 
plate. Thus, shot noise, arising in the galaxy survey, 
leads to a biased estimate of the amplitude. Perturb- 
ing the ISW template by T/ sw = Tisw + leads to an 
underestimate of the amplitude. Given shot noise vari- 
ance, a 2 /Ng a i, that scales with the number of galaxies, 
the expected amplitude takes the form, 



(A') = (A) 



1 



2 /N gal - 



(9) 



In the case of the ISW effect, the covariance ma- 
trix is dominated by the primary CMB. We estimate 
the covariance matrix using the best-fit angular power 
spectrum from WMAP5 convolved with the pixel win- 



dow function and instrument beam: 



Cij 



^CiWfBfPitcosOij), where Wi is the pixel win- 
dow function, and B\ is the instrument beam. Pixel 
noise is added according to the noise properties of the 
map considered. In the case of the ILC map, pixel noise 
is neglected because it is not characterized. 

A correction to the covariance must be added to ad- 
dress the fact that we estimate the mean of the dis- 
tribution from the sample itself. The covariance is 
C[a — ((xi — x)(xj — x)) where x is the mean pixel value 



l 



becomes 
l 



N pix Yi x i- Expanding, this 

AT Ya( X i X a) ~ TV Y a ( XaX j) 



within the survey: I 

{XiXj 

7v2 Y a Yp( x a x p)' Although this correction is < 5% for 
the SDSS survey area, neglecting it leads to an overesti- 
mate in the significance of the ISW signal. 

Finally, the computed covariance matrix is ill- 
conditioned, leading to an unstable inversion. This is 
generally the case for a map covering only a fraction of 
the sky. We find a usable pseudo-inverse by computing 
the singular value decomposition of the matrix and zero- 
ing the noisiest modes. The eigenvalue limit was chosen 



such that a stable inversion was produced for all CMB 
maps we considered. 

4.2. Results 

The calculation is carried out over the intersection of 
the galaxy survey footprint and CMB foreground mask, 
in total 9298 pixels. 

We tested ISW maps generated from four redshift cuts 
listed in Table [I] fit to the MCMC CMB map. The error 
was determined with Eq. We find that the signal 

degrades in the wider redshiit ranges in samples III and 
IV. The maximum signal was found in sample II, which 
includes ~ 2/3 of the galaxies about the median redshift, 
spanning 0.48 < z < 0.58. Including all galaxies reduces 
the significance by 0.6<j. We attribute this to contami- 
nation from a large-scale gradient in the galaxy density 
arising from a systematic shift in the redshift distribu- 
tion discussed in §2.1. This error is not accounted for 
in our analysis and we limit ourselves to the sample II 
redshift range for our results. The differences between 
the ranges that we test are comparable to the size of the 
photometric redshift error, so adjusting the cut should 
not significantly bring in new structures affecting the in- 
trinsic signal. 

The power spectra of the ISW maps are shown in Fig. 
[4] The maps are multiplied by the best-fit A value from 
Table [l] This normalization properly scales the spectra, 
aligning the peak at multipole I = 4. This scaling prop- 
erty with A suggests that narrowing the redshift range 
acts as a smoothing kernel, washing out fluctuations in 
the potential. We further investigate this effect with sim- 
ulations in §4.3. 

We calculated the best-fit ISW amplitude for the Q, 
V and W maps, as well as for the various foreground 
reduced maps using sample II. The results are listed in 
Table |2j and are all consistent with a 2a signal. 

There is little evidence for foreground contamination. 
Although there is a 0.2a difference between the V and W 
maps, the amplitudes found from the foreground-reduced 
maps agree with the raw co-added maps. This suggests 
that the slight color dependence is not due to foreground 
contamination, but may arise from other differences such 
as the map beams and noise properties. Furthermore, the 
ILC and MCMC maps agree with each other, and match 
the V filter. This finding is consiste nt with previous ISW 
studies, including |Ho et al.| (|2008| ) who place an upper 
limit of 0.3a on tne effect of foregi 



eground contamination 



on the ISW detection significance. 
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Li et al. (2009) found a large-scale correlation in 
WMAP data between the number of observation passes 
made over a region of sky and the derived temperature 
which could affect the ISW measurement. To check for 
this, we computed the template fit between the number 
of observations of each pixel and the temperature of the 
Q map at n S id e = 64 resolution. Over the SDSS area, we 
find no correlation on the mean pixel temperature with 
A = (6 ± 40) x 10~ 6 fjK/N oba . It is still possible that 
there is an effect from the disparity in observation num- 
ber between th e plus and minus antennae as reported by 
Li et al. (2009) which will require further investigation. 

4.3. Interpretation 

The constructed potential is affected by the survey 
window function, photometric redshift errors and shot 
noise. We account for these effects with the aid of simu- 
lations. We generate mock LRG catalogs from Gaussian 
large-scale structure simulations extending to z = 2 with 
40-Mpc cells and model the SDSS LRG sample as a slab 
cut out from this volume. We used a fast fourier trans- 
form to compute the potential in this rectangular ge- 
ometry, enabling the results of hundreds of simulations 
to be averaged together. The template fit in this case 
measures the contribution of the survey slab's ISW sig- 
nal to the full ISW signal integrated from z = to 2. 
We model redshift errors by convolving the density field 
with a Gaussian kernel with <j z = 0.04 along the radial 
direction and introduce shot noise by sampling the field 
according to the survey selection function. 

Fig. [5] shows the effect of shot noise in the simulations 
and on the LRG measurement. We under-sampled the 
LRG catalog to create ISW maps with increased shot 
noise and fit the trend with Eq. d9k. Although the am- 
plitude measured from LRGs is higher than from sim- 
ulations, the trends are similar. Shot noise biases the 
measurement low; in Fig. [5j the asymptotic limit to A 
from the LRGs is 2.8, to be compared with the measured 
result from the full map, 2.51, in Table |2j 

The redshift error convolution and survey window 
function act to bias A high. This can be understood 
due to the smoothing effect redshift errors have on the 
reconstructed potential. Fig. [6] shows the results of the 
simulations using the four tested redshift cuts. The mea- 
sured trend in A with survey width does not agree with 
the simulation results, but we attribute the deviation in 
the wider redshift samples to the large-scale gradient in 
the galaxy density in the tail of the redshift distribu- 
tion (discussed in §2.1), and so, we do not consider these 
points for the cosmological interpretation. 

After modeling both shot noise and redshift error ef- 
fects, we make a simple correction to our measurements 
by normalizing A by the results measured from simula- 
tions. In sample II, A is corrected by a factor of 1.2. The 
corrected values are listed in Tables [T] and |2j Our mea- 
surement is a factor of ~ 2.1 above the ACDM simulation 
results, within the 1-cr range of cosmic variance. This 
result is in concordance with the higher - than-expected 
ISW amplitu de found by Ho et al. (2008), Giannantonio 



et al. 



corre. 



([2008]) and others, and is consistent with our own 
at ion function measurements presented in §5. 

5. CROSS-CORRELATION AMPLITUDE 

5.1. Galaxy - CMB temperature correlation 



TABLE 1 

Dependence on redshift cuts 





redshift 


N 


Amplitude 


Amplitude 






range 


/10 3 


uncorrected 


corrected 


a 


I 


0.49-0.56 


300 


3.12 ± 1.63 


2.2 


1.9 


II 


0.48-0.58 


400 


2.51 ± 1.25 


2.1 


2.0 


III 


0.47-0.60 


500 


1.85 ±0.97 


1.7 


1.9 


IV 


0.45-0.63 


600 


0.95 ±0.66 


1.0 


1.4 



TABLE 2 
Amplitude fits 





Amplitude 


Amplitude 




Map 


uncorrected 


corrected 


a 


Q Coadd 


2.40 ± 1.19 


2.0 


2.0 


V Coadd 


2.51 ± 1.19 


2.1 


2.1 


W Coadd 


2.13 ±1.17 


1.8 


1.8 


Q FG reduced 


2.33 ±1.19 


1.9 


2.0 


V FG reduced 


2.51 ± 1.19 


2.1 


2.1 


W FG reduced 


2.20 ±1.17 


1.8 


1.9 


MCMC 


2.51 ± 1.25 


2.1 


2.0 


ILC 


2.51 ± 1.25 


2.1 


2.0 
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Fig. 4. — Plotted is a comparison of the power spectra of ISW 
maps constructed with different redshift samples, see Table [Tl The 
maps were multiplied by the best-fit amplitude, A. This factor 
properly normalizes the maps, which otherwise show a variation in 
amplitude. The signal primarily arises from low multipoles, / < 10. 

The standard measure of the ISW effect is the galaxy 
overdensity-CMB temperature cross-correlation. The 
measurement has be en made on the SPSS LRG dataset 
at th e 2 — 2.5a level ( |Ho et aL]|2QQ8 Giannantonio et al. 
[2008] ). 

The expect ed angular power spec trum in the flat-sky 
limit, see e.g. Afshordi et al. (2004), is an integral over 
the line-of-sight distance r, 



cr- 



: ^CMB 



3HQQ m bg 



J dr r 2 n(r) 



(l + 1/2)2 
d(l + z)D 1 (z) 
dr 



1 + 1/2 



(10) 



d fajiy is the galaxy selection function 

J dN 
dzdV 



where, n(r 
normalized by J r 

P (z,k = l+1 J 2 ^j is the matter power spectrum at red 



dr, Di is the growth factor and 



shift z. The correlation function in real space is, 
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Fig. 5. — Shot noise in the galaxy survey leads to an underes- 
timation of the amplitude. Plotted is A measured from a galaxy 
field with the given projected number density. The expected value 
of A is 1. The square markers are mean measurements of the ISW 
amplitude made on the under-sampled LRG catalog, which, when 
fully sampled, contains 53 galaxies/square degree. The error bars 
represent the sampling error on a single realization. The curves are 
fits to Eq. |9|. The solid curves represent mean results from mock 
LRG catalogs drawn from Gaussian simulations, with photometric 
redshift errors modeled by a Gaussian kernel with dispersion <r z . 
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Fig. 6. — We simulate the effect of redshift cuts on the de- 
rived ISW amplitude A (solid lines). The square markers show our 
measurement on LRGs, as in Table ^ We find that introducing 
redshift errors boost the measured amplitude by a factor of 2. The 
effect of shot noise is a downward bias in A, making the expected 
value less than 1 when cr z = 0. The dashed curve is the simulation 
result multiplied by 2.1, demonstrating that our measurement is 
a factor of 2 larger than expected in ACDM, although the mea- 
surement deviates from the simulation trend in the wider redshift 
ranges due to systematic redshift errors (see text). The error bars 
on the simulations represent sampling error on a single realization. 



where Pi are the Legendre polynomials. 

We measure the goodness-of-fit against a fiducial 
model scaled to the data. A one parameter fit is suf- 
ficient because variations in cosmology affect the ISW 
amplit ude, but have little effect on the shape of the spec- 
trum (Padmanabhan et al. 2005). We apply the s ame 
templ a te fit estimator (Eq. ^) in the fashion of [Ho 



et 



al. (2008): given an observed data vector, X, and 
moHel template, Y, we find the best fit amplitude, A, for 
the model X = XY. The maximum likelihood estimate 
of the scale amplitude A is A = yc^Y •> with variance 
a 2 = (YC^Y)- 1 , where C = (XiXj) is the covari- 
ance matrix. In the case of the angular power spectrum, 
C = (C/C//). Given a fixed cosmological model, the free 
scaling is a constraint on the product of the galaxy bias 



and normalizatio n of the power spect rum, b p cr|. 

We use SpICE ( Szapudi et al.|20 01) to estimate the C\ 
spectrum for multipoles / = 4 — 192, the limits of which 
are set by the resolution of the map and the extent of the 
survey. The measured angular power spectrum is binned 
into 15 logarithmic band powers. 

We consider only CMB variance in the error analysis. 
Neglecting the cosmic variance in the gala xy field leads 
to an underestimate of the errors by 10% (Cabre et al. 
2007). We estimate the covariance matrix in a Monte- 
Carlo procedure with 2000 realizations of the CMB gen- 
erated according to the best fit WMAP5 angular power 
spectrum. We consider here the correlation with the 
MCMC CMB map, which we model using an instrumen- 
tal beam with 1° full- width half-max. 

The cross-correlation results for the MCMC map are 
plotted in Fig. [7] The best fit amplitude is b g a\ — 
2. 82 ±1.35, a 2.1<j measurement. This is consistent with 
previous results using the LRG datasetplo et al. 2008 



|Giannantonio et al. 2008[ ). We find that the signal is 
a factor of 1.9 greater than the ACDM prediction with 
b g =2.2 and a 8 0.817. 

5.2. ISW-cleaned CMB map 

We now ask whether the constructed ISW map 
contains the signal found in the galaxy-CMB cross- 
correlation. We produce a cleaned CMB map by sub- 
tracting the constructed ISW template. The template 
is scaled by the best fit amplitude from Table [2] The 
cross-correlation of the cleaned CMB map with the LRG 
map is shown in Fig. [8j We find that the subtraction 
effectively nulls the correlation signal. This agreement 
strengthens the claim that the correlation is of the form 
expected from the ISW effect. 

5.3. ISW - CMB temperature correlation 

The ISW power spectrum gives a further measure with 
which to characterize our ISW construction. The cross- 
correlation between the reconstructed ISW map and the 
CMB temperature measures the ISW auto-power spec- 
trum, and, by construction it is boosted by one power 
of the galaxy bias: b Q Cj T . In the fla t-sky limit the ex- 



pected spectrum, e.g. Cooray (2002a), is, 



TT 



CMB 



(I + 1/2) 4 



I 



dr r 



d{l + z)D 1 {z) 
dr 



1 + 1/2 



(12) 



The ISW power is dominated by structures at low red- 
shift which are not captured by our map. To account 
for this, we limit the integral in Eq. (12) from z = 
0.48 — 0.58, the approximate redshift range of the survey. 
We carry out a fit to this model, again using Monte-Carlo 
realizations of the CMB to estimate the covariance. The 
measurement is plotted in Fig. [7| The best-fit amplitude 
is a factor of 1.74 ± 1.46 above ACDM. Though the de- 
tection is marginal, it is consistent with the amplitudes 
measured above. 

6. SUPERVOIDS AND SUPERCLUSTERS 

Recently we detected the imprints of voids and clus- 
ters on the CMB and attributed the finding to the linear 
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Fig. 7. — The ISW signal may be measured with the CMB temperature - galaxy cross-correlation function, plotted in panel a, top. 
The left frame shows the cross power spectrum in harmonic space, which corresponds to the real-space correlation function on the right. 
Over-plotted are ACDM template ISW models: the dashed blue line is the template scaled by the best fit amplitude and the solid red line 
is the ACDM prediction with erg = 0.82 and b = 2.2. Panel b shows the correlation of the constructed ISW map with the CMB, which 
measures the ISW auto-power spectrum. 

CMB are plotted in Fig. [9j We re-examine the signifi- 
cance of the correlation with the template-fitting analysis 
used above, using a template constructed with the mean 
profiles of the voids and clusters. A compensated model 
was chosen to fit the profiles to ensure that the mean 
of the map is zero. For simplicity we use a "Mexican 
hat" Laplacian of a Gaussian for the functional, fitting 
for the amplitude and width. The template fit confirms 
the 4a significance of the measurement found previously: 
we find A = 1.2 ± 0.25. This high signal is not surprising 
because the void and cluster profiles were measured on 
the CMB, although it does confirm the 4a peculiarity of 
these sites on the CMB. 

If this correlation is due to the linear ISW effect, we 
would expect the signal to be contained in the con- 
structed ISW map presented here. However, we find 
that the mean temperature of the clusters and voids on 
the ISW map is not significant. The temperature dif- 
ference between clusters and voids is T cluster — T volds = 
0.08 ± 0.1/iif , and stacking the clusters and voids sepa- 
rately does give weak hot and cold spots. We measure 
the temperature in a 4° compensated aperture. The er- 
ror bar was measured through a Monte-Carlo procedure 
using random locations within the survey. This result 
contradicts the suggestion that the signal is due to the 
linear ISW effect. 

7. CONCLUSIONS 
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Fig. 8. — Subtracting the constructed ISW anisotropy from the 
CMB map nulls the correlation. The solid curve shows the original 
galaxy-CMB cross correlation. After subtracting the ISW signal, 
it is reduced to the dashed curve. The measurement error is the 
same as in Fig. [7] 



ISW effect ( |Granett et aTl|2Q08| >. The same SDSS LRG 
dataset at z = 0.5 discussed here was used in that work, 
although with minor differences in the selection criteria 
and photometric redshifts. Thus, with the ISW signal 
reconstruction, we can further investigate the result and 
place limits on the role of linear ISW in the measurement. 
We also reassess the significance of the measurement with 

a template fit on the CMB. 

The 50 voids and 50 clusters presented by |Granett et al.| 
(2008) are strongly correlated with the CMB. The clus- 
ters, on average, fall on hot spots on the CMB, while the 
voids are cold. The mean temperature profiles on the 



We implement a template fit measurement of the ISW 
signature imprinted on the CMB from structures traced 
by SDSS luminous red galaxies at z ~ 0.5 with photo- 
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Fig. 9. — Left: the measured void and cluster temperature pro- 
files on the CMB from jGranett et al.||2008| >. Right: we use the 
mean profiles to generate a map ot the voids and clusters on the 
sky. Voids are blue and marked by a solid circle; clusters are red 
and marked by a dashed circle. We assess the significance of these 
temperature fluctuations with a template fit, confirming the origi- 
nal finding. 

metric redshifts. The approach uses the density informa- 
tion in the galaxy survey to remove the effects of local 
cosmic variance from the ISW measurement. The de- 
tection significance is 2<r, which is consistent with the 
cross-correlation statistic. 

Our measured amplitude confirms the cross-correlation 
finding, that the signal is la higher than the ACDM pre- 
diction. The evidence suggests that this am plitude arises 
from variance on the CMB ( |Ho et aL||2008| ), or alterna- 
tive ISW physics, although we cannot exclude cosmic 
variance in the galaxy sample. 

In principle, the template fit simplifies the error anal- 
ysis and cosmological interpretation, because a Monte- 
Carlo procedure is not required to assess the role of cos- 
mic variance in the galaxy sample, and multiple datasets 
can be combined naturally. On the other hand, the 
measurement is influenced by a number of systematic 
effects that do hinder the way towards precision cosmol- 
ogy. First, contamination by foreground emission must 
be subtracted from the CMB and the accuracy of the 
subtraction ultimately limits the ISW detection. This is 
also true of cross-correlation function studies. Second, 
shot noise in the galaxy survey biases the estimator, but 
in a largely correctable fashion. In the case of the LRG 
survey, we estimate that the effect is < 10%. Third, the 
photometric redshift errors severely degrade the recon- 
structed density and potential 3D maps, leading to a bi- 
ased estimate of the template amplitude. In simulations 



of the LRG catalog, the template signal was reduced by 
a factor of 2. We expect that this bias can be mitigated 
by using cosmological density reconstruction techniques 
that take the red shift errors as a Bayesian prior (e.g. 
Enfilin et al.||2QQ8] ), and the effect can be corrected with 
simulations. 

The ISW template fit provides a unique constraint on 
cosmological parameters. In particular, because the ISW 
reconstruction is made from a biased tracer, the template 
fit constrains the galaxy bias independently of the the 
power spectrum amplitude, <r 8 . Precision results from 
the ISW signal will become feasible with new all- s ky sur- 
vey projects, such as Pan-STARRS flKaiser||2QQ4l ) . 

In this context, we further investigated our measure- 
ment of imprints of ~ 100 Mpc supervoid s and super- 
clusters on the CMB ( |Granett et al.| [2008). We confirm 
that the signal is present in WMAP data at 4<j confi- 
dence with a template fit. However, we find that the 
signal has much lower amplitude in the linear ISW map 
constructed from LRGs than in the WMAP data. 
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